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We discuss the most interesting approaches to derivation of the Bekenstcin-Hawking 
entropy formula from a statistical theory. 

in' 

O ' 1. Introduction 
O 

' Black hole is one of new physical phenomena predicted by General Relativity (GR). Though 

the existence of black holes is still hypothetical (see, for example, [1]), they send a challenge to our 
imagination and rise a lot of interesting questions. The space-time having a black hole in it, first, 
\ has a singularity, and second, has a horizon preventing an external observer from seing it. The 
singularity in GR is radically different from field theory singularities because it is a property not 
of some field but of the space-time itself. The topology of space-time is changed when it acquires a 
^ \ black hole. 

I The most convincing viewpoint is that arising of singularities demonstrate the incompleteness of 

classical gravitational theory and the need to replace it by the quantum theory. So, if one believes 
in GR as a right theory of gravitation then he can assert that the singularity problem will be solved 
\ only after discovering of the quantum gravity theory. But. . . might quantum gravity be the End 

■ of Physics? If the problem of black hole singularity would be solved then probably the problem of 
O . cosmological singularity will be solved also. And so, we will have an answer to the question how, 

when and from what kind of material our Universe has been created. Then all physicists together 
should migrate, maybe, to biology. 

But sensible people suggest that the End of Physics is still far from us, and so, we are still far 
from quantum gravity. Then what we can expect now? Maybe, one should try to guess at least one 
^ ' formula of the future mysterious theory, like Johann Jakob Balmer [2] , Swiss teacher of mathematics 

■ at a secondary school for girls in Basel, has done in 1885. Or, if one believes that such a formula has 
been already discovered by Bekenstein [3, 4, 5], one should hope to derive it from some intermediate 
between classical and quantum theories model, like Niels Bohr has done in 1913 with his model of 
hydrogen atom. 

Thorough study of black hole solutions and their properties such as test particle behavior has 
led to formulation of "black hole thermodynamics" [6], i.e. four laws which are very similar to four 
laws of the thermodynamics: 

0. Surface gravity k is constant over the horizon. 

1. Energy is conserved. 

5M = ^ + nSJ + ^5Q 

2. Matter entropy plus sum of areas of all horizons can not decrease. 

-r {Smatter + Sbh) > 

at 

3. It is impossible to reduce surface gravity to zero by a finite sequence of operations. 

This discovery was motivated by an observation that black holes revealed irreversibility: a 
particle falling into the black hole was unable to be back. This particle enlarges the black hole mass 
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and also the gravitational radius and the horizon area. It is impossible to decrease these quantities 
in classical physics [7, 8]. 

The entropy is an example of non-decreasing variable in the standard thermodynamics. Beken- 
stein [3, 4, 5] attributed entropy to black holes by suggestion that it is proportional to the horizon 
area. The constant of proportionality was determined by Bekenstein and Hawking in articles [5, 9, 
10] . So, the following formula has appeared 

Sbh = A/ AG = A/Uli > (1) 



where Sbh is the Bekenstein-Hawking entropy and Ipi = ^JhG/c? is the Planck length. 

Another characteristic variable of thermodynamics is the temperature. Black holes also have 
analogous variable. This is the surface gravity n which (in case of a Killing horizon) is defined by 
the following equation 

X^^cXb = KXb, (2) 

where is a null Killing vector tangential to the horizon. Practically this temperature is rather 
low for black holes which could be formed as a result of a star collapse. If Mq is mass of the Sun 
then in standard Kelvin units 

Teh = K/27T = -4— M"^ « IQ-^Mq/M. (3) 

OTTKB 

Nevertheless Hawking suggested that taking a black hole together with radiation of the same 
temperature and putting them into a thermostat should provide a system in thermodynamical 
equilibrium. Therefore the black hole absorbing radiation have also to emit. Under some conditions 
Hawking has predicted [10] that a process of emitting particles takes place for any black hole. Here 
the gravitational field of black hole is treated as a classical background whereas the field of emitted 
particles is quantized. 

The Hawking radiation becomes relatively stronger if the black hole becomes smaller. But the 
back-reaction of emitted particles on the gravitational field was neglected. Hawking had enough 
courage to say that a black hole could even explode or evaporate. This suggestion has led to a new 
puzzle: if a black hole disappears then what can be said about the information "eaten" by this hole? 
What if unitarity or even quantum mechanics are violated by this process? This is an excited and 
open question (sometimes called "black hole information puzzle"), but we will not discuss it here 
(see [11]). 

The subject of this talk is a review of statistical derivations of the black hole entropy formula. 
History of physics shows that all thermodynamical laws later have been derived from a more fun- 
damental theory — statistical mechanics. It was shown that the concept of temperature could be 
explained as the average kinetic energy of a micro-particle, the concept of entropy — as logarithm of 
the number of states corresponding to the same macroscopic thermodynamical macro-state. Look- 
ing at the black hole from this viewpoint we are to relate its entropy with the number of states of 
some quantum theory. If we consider Bekenstein-Hawking formula in some sense analogous to the 
Balmer formula for hydrogen atom then it is required to construct a phenomenological model of 
quantized black hole which would lead to Eq. (1). 

It is natural that there are different approaches to this problem. Which one is true or maybe 
all of them are misleading, nobody knows now. Below we demonstrate in short the most popular 
methods and their achievements. Personal sympathies of the author can be seen from proportions 
of text devoted to one or another approach. 
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2. Conformal field theory 

2.1. Thermodynamics of almost extremal black holes 

There is an important analogy described by Fursaev [12]: thermodynamical properties of almost 
extremal black hole correspond to those of a gas of massless scalar particles in flat 1+1-space. We 
start with explanation of this idea below. 

A charged black hole is described by Reissner-Nordstrom solutions 

ds"^ = -Bdt^ + ^ + r'^dn'^, (4) 
B 

where dO? is the metric on a unit sphere, 

B = —^{r — r-){r — r+), r-t = m ± \/ m? — q^, (5) 

q = QVG is related to the electric charge Q of the black hole, and m = MG to its mass M. The 
radius of the horizon is r+. The Hawking temperature of this black hole is 

TH = -^^fi^^^, (6) 
27rr+ 

and the Bekenstein-Hawking entropy is Sbh = 'Kr\/G. 

This solution has an interesting property: the Hawking temperature vanishes when m ^ q. A 
solution with m = q is called an extremal black hole. There are no physical processes allowing 
to convert a black hole with m > q into an extremal one. Nevertheless the extremal black holes 
acquired a lot of theoretical attention. 

Let us consider now black holes which are "almost extremal" 

m = q + E, E'€.q. (7) 
Thermodynamical relations for these objects are 

E = m-q = X''T^H, S = Sbh - = ^T„, (8) 

where A = (27r2g3)V2, 

Now consider a gas of massless non-interacting scalar fields (f)k, k = 1, ... ,c on an interval [0, 6]. 
The dynamical equations and boundary conditions are as follows 

{dl-dl)Mt,x) = o, Mt,o) = Mt,b) = o. (9) 

Let this system be in a state of thermal equilibrium at some temperature T. Then it is a one- 
dimensional analog of an ideal gas of photons in a box. The free energy of this gas 

F{T, L) = cTY^ In (l - e-'^"/^) , (10) 

n 

in the thermodynamical limit, T6 S> 1, can be estimated as 

F{T,h)^-'^bT\ (11) 

so, the energy £^(r, h) and the entropy S{T, b) of this system are the following 

E{T,h)c,^bT\ SiT,b)c^^bT. (12) 

By comparing (12) with (8) it is easy to see that thermodynamical properties of an almost extremal 
black hole are identical to properties of an ideal massless gas in a flat two-dimensional space-time. 
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2.2. BTZ black holes 



The first statistical derivation of the black hole entropy formula based on the conformal field 
theory (CFT) was suggested by Strominger [13]. We will follow him below in this subsection. 
Consider 3-dimensional GR with a negative cosmological constant A = — Z"^: 

One of solutions for the corresponding gravitational equations is anti-de Sitter AdS^ space: 

„2 _ / 1 , \ Ji2 , / 1 , ^ \ . 2 , „2 ...2 



ds' = - (^1 + -J dt' + (^1 + dr' + V, (14) 

where 99 has period 27r. There is another exact solution discovered by Bahados, Teitelboim and 
Zanelli [14] and later called the BTZ black hole (here we consider only non-rotating case): 

Brown and Henneaux [15] studied a group of the asymptotic transformations, i.e. the transforma- 
tions preserving the property of a metric to be "asymptotically anti-de Sitter" . The BTZ geometry 
is asymptotically AdSz and fulfils the following boundary conditions first given in paper [15]: 



gtt 


^2 

= 1^ + ^(1 


gt(i> 


= c?(i), 


9tr 


= ^(^)' 




f 1 


9rr 




9r4> 




9H 


= r2 + 0(l) . 



(16) 



Strominger has shown [13] that at large r the diffeomorphism vector fields 5x^ = C^{x) preserving 
this asymptotic structure (16) are as follows: 

= £(T+ + T-) + ^{dlT+ + dlT-) + 0(r-^), (17) 

C^^ = (r+ - T-) - ^{dlT- - dlT~) + 0(r-^), (18) 

C = -^r(9+r+ + 9+T+) + O(r-i), (19) 

where d± = £dt ± d^p, T"*" and T~ are functions of one variable, t/l + ip and t/l — ip correspondingly. 
The algebra of these vector fields, where the commutator, or a Lie bracket, [Ci,C2] of two vector 
fields, and ^2 > is a vector field with components 1^3 = C,idi,C,2 ~ 1 preserves the asymptotics. 
The commutator of vector fields Ci and C2 with the above asymptotic behavior is a vector field Cs 
which has the same asymptotic behavior, and = T^d±T^ — T^d±T^. So, generators of these 
diffeomorphisms should form a closed algebra. Due to the periodicity over variable (p and the 



57 



specific dependence of functions on t these functions also are periodic over time coordinate with 
period 27r£. 

Let Cnj Cn be defined by relations (17)-(19) with 

r+ = |e^"(*/^+^), r- = o, r+ = o, f- =^-e'<^'^-^\ (20) 

The algebra of these vector fields is as follows 

[Cn, Cm] = (n - m)Cn+m, [Cn, Cm] = (n - m)Cn+m, [Cn, Cm] = 0. (21) 

and therefore the asymptotic symmetry group (17)-(19) generated by two copies of the Virasoro 
algebra is the conformal group. So a consistent quantum gravity theory on AdS^ should be a 
conformal field theory. It is well known that due to conformal anomaly the Virasoro algebra has a 
central extension. Surprisingly Brown and Henneaux have derived this result even classically from 
surface terms appearing in Poisson brackets of the gravity Hamiltonians. 

In the Hamiltonian formalism of GR any infinitesimal diffeomorphism 5x^ = C,^ (x) is generated 
by a Hamiltonian of the following form 

^[C]= / C^W^ciS + JiC]. (22) 

This expression is a sum of a volume integral of a linear combination of constraints Ti^ = with 
the diffeomorphism parameters as their coefficients, and a surface integral depending on the same 
parameters and on canonical variables, gij, and 'Kij. The surface term J[Q is introduced to ensure 
the canonical form for variations of i?[C], 

SH[C]= [ {A'^gij + B'^diT^,), (23) 

(without any surface terms) according to the ideology of Regge and Teitelboim [16]. Both the 
very existence of J[C] (so, the very existence of i?[C]) and the form of it depend on the boundary 
conditions. When the constraint equations are fulfilled, H[(] reduce to J[C]- For the BTZ black 
hole solution the value of H coincides with the black hole mass M. 

Let us define generators L„ and L„ by relations L„ = i?[Cn], -^n = -ff[Cn]- The Poisson brackets 
of these generators were calculated by Brown and Henneaux [15]: 

[Ln, Lm] = (n - m)Ln+m + J^i^ ~ ^)<^n+m,0, (24) 
[Ln, ^m] = {n- m)Ln+m + J^i'^^ ~ n)Sn+m,0, (25) 

[Ln,Lm] = 0, where the constant c is the central charge given by 

(26) 

For the BTZ black hole Bekenstein-Hawking entropy is 



where M = r\/{Sl'^G). This formula looks like the formula (12) for the CFT entropy given in 
previous subsection, where extension I in Eq. (27) is analogous to parameter h in (12). 
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There is the following relation between the energy and the Virasoro generators 



H = j{H[Co] + H[Co]) = \{Lo + Lo). (28) 

Let up to constraints Lq = h, Lq = h. If the energy of the system coincides with the mass of the 
black hole, then, due to h h symmetry, h = h = Ml/2. 

The main suggestion by Strominger is the existence of a quantum theory of gravity on AdS^ 
and the coincidence of its results with the results of classical theory given above. 

In the semi-classical regime the spectrum of a black hole can be treated as continuous. The 
classical black hole is a highly degenerate state /i + ^ » c » 1. The degeneracy D, D of operators 
Ln, Ln , can be found by using the Cardy formula [17] 




InD ~ 27r^/— = 27r;/-^. (29) 

The total degeneracy is 



ch W^M 

InD + lnD = 21nD ~ 47rW— = 27rW— — , (30) 

V 6 V 2G 

which is exactly the Bekenstein-Hawking entropy of the BTZ black hole with mass M. 

An obvious drawback of this derivation is certainly that it is limited to the 2 + 1-gravity case. 
Though the idea of transformation of gauge degrees of freedom living on the boundary into dynam- 
ical variables due to boundary conditions is attractive, but in the 2 + 1-dimensional space-time it is 
irrelevant what part of the boundary (outer or inner) is used while in 3 -|- 1 case this is already not 
so. Therefore boundary conditions should be specified on the black hole horizon and not on spatial 
infinity. 



2.3. Carlip I 

Carlip [18] invented a new method for statistical calculation of the black hole entropy in order to 
apply the ideas by Strominger [13] discussed in the previous subsection to a general n-dimensional 
space-time, without the restriction to 2 -|- 1-dimensional case. 

The metric of n-dimensional space-time may be taken as follows 

ds'^ = -N'^dt^ + f{dr + N'-dtf + a^pidx'^ + N'^dt){dx^ + N>^dt), 

where a,f3... correspond to coordinates on a sphere r = const, t = const and run over values 
1,2,... ,n — 2. The function tends to zero on the horizon of a black hole r = r+ in such a way 
that 

N'^ = h{x''){r-r+) + 0{r-r+f. (31) 

Further, the Hamiltonian looks like a linear combination of constraints {'Ht,'Ha} plus, probably, 
some surface integrals 

Hoii] = I d^-^xen^, m = Ho[i] +<f ... , (32) 

where the parameters of deformations of a constant time hypersurface (Lagrange multipliers for 
constraints) are components of the decomposition for an infinitesimal space-time diffeomorphism 
over the basis {{l/N,-N''/N),d/dx''}: 

= Ne, e = r + iv"e*- (33) 
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If one calculate a variation over gah and tt"'' of that part of a Hamiltonian containing only constraints 
(without the boundary terms), and consider as not depending of canonical variables, one get 



SHo = S 



(34) 



where ra" is a unit normal to the boundary at t = const, K^b is an extrinsic curvature tensor of a 
constant time hypersurface, tt"^ = —j\fo{K°^ — g°^K) are the momenta conjugate to the spatial 
metric gab- The surface integral is taken over a boundary which should include both the horizon of 
a black hole, and the spatial infinity. No boundary conditions arc applied here. Up to notations 
the formula coincides with the similar formula from Reggc and Tcitclboim work [16]. 

Then it is necessary to put boundary conditions for the canonical variables and for the parameters 
of deformations. They are 

1. to set in a phase space a domain "close" to a black hole solution; 

2. to provide preservation of this domain under evolution determined by the parameters of de- 
formations. 

Carlip [18] starts with the following boundary conditions at the horizon: 

/ = ^iV-' + 0(l), W = 0{N\ a„;3 = 0(l), iV" = 0(l), 

(dt - N'-dr)g^, = 0{N)g^,, VaNp + V^jAT^ = 0{N), drN = 0{1/N), 
Krr = 0{1/N^), K^r = 0{1/N), K^p = 0(1), 

f = o(iv2), = o(iv), r = o(i). 

What about a part of the boundary arranged at spatial infinity we get rid of the corresponding 
contribution to the Hamiltonian variation by considering parameters ^ as rapidly decreasing when 
r — >■ oo. 

Despite of the accepted boundary conditions the variation of Hamiltonian Hq[^] still contains 
some surface contribution and to get rid of it is suggested [18] to add a surface integral of the 
following form 

J[i\ = ^ jf_ d"-'2;{n«V„|* + + [uatK^] } . (35) 

Later the extra restrictions K„- = = K^p are imposed on the exterior curvature tensor. As a result 
the last term in the variation formula (34) becomes exactly zero and, in fact, the term written in 
square brackets of Eq. (35) is not required by the Regge-Teitelboim ideology, just opposite, it spoils 
the "differentiability" condition (in [18] this is masked by an unjustified restriction SKrr/Krr = 
0{N)). So, the variation of the improved Hamiltonian L[^] = i?o[?] + <^[^] becomes 

6L[^] = bulk terms I dT'^x 5rfdr£,*^. (36) 
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Then it is argued [18] that, first, we can calculate a Poisson bracket for generators of two various 
deformations according to the following formula 

{L[i2],m]=d^^L[i,], (37) 

and second, that this equality remain valid after reduction, that is, after putting constraints to 
zero "in strong sense". It means, that the formula (37) remains valid for surface integrals taken 
separately: 

boundary terms o/|l[^2], ~ boundary terms of S^^Hq[^i] + S^^J[^i]. (38) 

In the same time the Poisson bracket (37) should have the following form 

d^^Liii] = = L [{6,6}sd] (39) 

where if [^1,^2] is a possible central term in the algebra. Here {^i,^2}sD is the Lie bracket for the 
algebra of surface deformations, given by 

Eq. (39) will be used to compute the central charge. 

As a result, the following expression was obtained for the surface contribution to the Poisson 

bracket {L[e2], 

" 8^ ^"'^^ {72 [dr{ft2)drii - dr{ftl)drii] + jdr " S^J\\ | . (41) 

In further it is used to realize the Virasoro algebra which central charge. 

Then Carlip [18] reduce the general formalism to the case of an axially symmetric black hole, 
with an adapted angular coordinate (p such that dtpg^i, = and consider a particular subalgebra of 
the surface deformation algebra with the following properties: 

1. The surface deformations are restricted to the r-t plane. The r-t plane has the central role 
in determining the entropy. 

2. The difFeomorphism parameter ^* = ^*/iV "lives on the horizon," in the sense that near r = r+ 
it depends on t and r only in the combination t— r*, where fdr = Ndr^, in the time-slicing such 
that Nr = 0. For the Kerr black hole, t — r^, is essentially the standard Eddington-Finkelstein 
retarded time, up to corrections of order r — r+. 

3. The lapse function N is fixed at the horizon. The horizon is physically defined by the condition 
N = 0, while our boundary term (35) is written at r = r+; this condition ensures that the 
boundary remains at the horizon. 

Condition 1 imply that the diffeomorphism parameter ^'f' has the form 

N't' ^ 

Condition 2 requires that 

dre = -j;^dte (43) 
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at r = r+, allowing us to write radial derivatives at the horizon in terms of time derivatives. To 
impose condition 3, we can examine difFeomorphisms of 5** = — l/AT^. With initial coordinates 
chosen so that Nr = 

V* = ^ = ^{9t- N%) e + ■ (44) 

This structure suggests that the difFeomorphisms are separated into left-moving modes ^* with 
dt^* = fiS^^*, and right- moving modes ^* with 5^^* = —^d^^*, where Q = —N'^{r+) is the angular 
velocity of the horizon. Then 

e = —f^dte, e = o. (45) 

Now it is possible to write the left-moving modes at the horizon as 

C*=^exp|^(t-r. + n-V)}, (46) 

where T is an arbitrary period. The normalization has been fixed by the requirement that the 
surface deformation algebra (40) reproduce the Diff algebra 

{im, LYsB = - ^)Cm+n- (47) 

Substituting the modes into the boundary term (41), we obtain 

A j3 

(5|^L[|„] = hulk terms + '^:^f^^^^rn+n, (48) 

where A is the area of the boundary aX r = rj^. We can now use a trick of Brown and Henneaux 
to evaluate the central term K[^rn,in]- When evaluated on shell, the Hamiltonian and momentum 
constraints vanish, so Hq[^] = 0. Equation (39) thus reduces to a collection of boundary terms. 



in 6jn+n = J[{^m, Cn}sD] +K[^,n, Cn] 



SttGT 

= i{n - m)J[im+n] + K[U, in]- (49) 

It is easily checked that 



on shell. Hence 



A (3 

K[^m,^n] = -^—^-inir? - —^)5m+n, (51) 



the correct form for the central term of a Virasoro algebra with central charge 



(52) 



The results above imply that the quantum states that characterize a black hole horizon must 
transform under a Virasoro algebra with central charge (52). This is sufficient to permit the use of 
methods from CFT to count states. In particular, a CFT with a central charge c has a density of 
states p{Lq) that grows asymptotically as 



logp(Lo)~2W^, (53) 
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where c^s is an "effective central charge" . If the ground state is an eigenstate of Lq with eigenvalue 
zero — then Ces = c. Following Strominger, let us assume these conditions are satisfied in quantum 
gravity. Then from Eqs. (50) and (52), 

logp(Lo)~^, (54) 

giving the standard Bekenstein-Hawking entropy. In general, right-moving modes may make an 
additional contribution to the density of states, but the central charge for those modes vanishes, so 
Eq. (54) gives the full entropy. 

Carlip [18] himself has provided some critical comments to his derivation: 

• Like Strominger's derivation, this computation uses symmetry arguments to derive the behav- 
ior of any microscopic theory of black hole horizon states. This is both good and bad: good 
because it provides a universal explanation of black hole statistical mechanics, bad because it 
offers little further insight into quantum gravity. 

• It also seems plausible that the description of black hole entropy here is related to the picture 
of microscopic states as "would-be pure gauge" degrees of freedom that become dynamical 
at a boundary. The existence of a central charge is an indication that the algebra of surface 
deformations has become anomalous, and that invariance cannot be consistently imposed on 
all states. 

• The extremal black hole is typically characterized by a lapse function behaving as ~ 
(r — r+)^ near the horizon. Such a configuration satisfies the boundary conditions assumed 
here, but in contrast to the non-extremal result it gives J[^o] = 0, i.e. classically Lq = 0. 

This approach has also been criticized for some other drawbacks by several authors: [19, 20] 

• The improved Hamiltonian remains non-differentiable. 

• Standard Poisson bracket calculation gives zero central charge. 

• It is inapplicable to spherically symmetric case. The axial symmetry is required. 

2.4. Carlip II 

Being unsatisfied with the above mentioned contradictions in his derivation of the black hole 
entropy [18] Carlip proposed another approach [21] to calculation of the surface terms in Poisson 
brackets. It was based on a definition of Poisson brackets through the symplectic current and 
covariant phase space technique. 

First, let us consider a general diffeomorphism-invariant field theory in n space-time dimensions 
with a Lagrangian L [</>], where L is viewed as an n-form and (f) denotes an arbitrary collection of 
dynamical fields. The variation of L takes the form 

SL = E ■ S4> + d& , (55) 

where the field equations are given by E = and the symplectic potential Q[4>, dcp] is an (n — l)-form 
determined by the "surface terms" in the variation of L. The symplectic current (n — l)-form u is 
defined by 

io[<p,6i<p,S24>] = 6iQ[4>,62cl)] - 52@[<P,Si<l>], (56) 
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and its integral over a Cauchy surface C, 

J^[0,<5i0,<520] = [ a;[(/.,(5i(/.,(520] (57) 
Jc 

gives a presymplectic form on the space of solutions of the field equations. This space, in turn, 
can be identified with the usual phase space, and $7 becomes the standard presymplectic form of 
Hamiltonian mechanics. 

For any diffeomorphism generated by a smooth vector field one can define a conserved 
Noether current {n — l)-form J by 

3[^] = e[ct>,jr^4>]-^-i^, (58) 

where denotes the Lie derivative in the direction ^ and the dot • means contraction of a vector 
with the first index of a form. On shell, the Noether current is closed, and can be written in terms 
of an (n — 2)-form Q, the Noether charge, as 

3 = dQ . (59) 

Now consider a vector field and the corresponding generator of diffeomorphisms In the 

covariant phase space formalism, Hamilton's equations of motion become 

5H[C]=n[<P,5(l),C^(l)]. (60) 

It is easy to see that when (f) satisfies the equations of motion, 

u;[<j>,dcf>,C^<j>] = S3[^] - ■ &[<!>, dcl>]), (61) 

so 

m]= [ (Q[e]-^B), (62) 

JdC 

where the (n — l)-form B is defined by the requirement that 

5 [ ^B[<^]= / e-G[<^,,50]. (63) 

JdC JdC 

Given a choice of boundary conditions at dC, finding B is roughly equivalent to finding the appro- 
priate boundary terms for the Hamiltonian constraint in the standard ADM formalism of general 
relativity. 

The Lagrangian n-form of GR in n space-time dimensions is 

Lai...a„ = ^Q^Q^ai...an^J (6^) 

yielding a symplectic potential (n — l)-form 

eai...a„_i [g, Sg] = j^eba,...ar._, [g''V dgdeSg"') " Ve*^^'') • (65) 

The corresponding Noether charge, evaluated when the vacuum field equations hold, is 

Qai...a„-2[<7,^] = -Y^ebca,...an.,V''e- (66) 

In the absence of a boundary, the Poisson brackets of the generators form the standard 

"surface deformation algebra" , equivalent on shell to the algebra of diffeomorphisms. On a manifold 
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with boundary, however, the addition of boundary terms can alter the Poisson brackets, leading to 
a central extension of the surface deformation algebra. That is, the Poisson algebra may take the 
form 

{H[^i],H[C2]} = H[{Ci,C2}] + K[C,,U (67) 

where the central term K[^i, ^2] depends on the dynamical fields only through their (fixed) boundary 
values. 

Consider the Poisson brackets of the generators of diffeomorphisms in the covariant phase space 
formalism. Let and be two vector fields, and suppose the fields cf) solve the equations of motion 
(so, in particular, the "bulk" constraints are all zero). Denote by 5^ the variation corresponding to 
a diffeomorphism generated by ^. For the Noether current J[Ci], 

% J[6] = J[6] = 6 • dJ[6] + d{^2 ■ J[Ci]) = d [6 • m4>,c^,<l>] - 6 • L)] , (68) 

where the fact that dJ = on shell was used. Hence 



Jc 



IC 

= [ (6-e[<^,£5,</>]-6-e[</>,£j,<^]-6-(6-L))- (69) 

JdC 

Then Brown-Hcnncaux [15] trick is used. Since equation above is evaluated on shell, the "bulk" 
part of the generator on the l.h.s., which is a sum of constraints, vanishes. Hence l.h.s. can 

be interpreted as the variation S^^J[i,i], where J is the boundary term in the constraint. On the 
other hand, the Dirac bracket {-/[^i], ^^[^2]}* means precisely the change in J[^i] under a surface 
deformation generated by ^[^2]; that is, 

%^[6] = {J[a],J[6]}*. (70) 

Comparing Eq. (67), evaluated on shell, we see that 

i^[6,6] = %^[6]-^[{a,6}], (71) 

where % is given by (69). Given a suitable set of boundary conditions, this allows to determine 
the central term ^^[^1,^2]- 

In particular, for GR without matter sources, the Lagrangian L vanishes on shell, and r.h.s. of 
(69) can be estimated from (65). We get 

{J[U = T^l^ e,ea,..a„_, [e2'Vd(V'^e? - V^^f )- 

- ?JV,(V'^ei - ^=^2')] • (72) 

Consider an n-dimensional space-time M with boundary dM, such that a neighborhood of dM 
admits a Killing vector x"' that satisfies = 9abX"'X^ = at dM. In practice, it is useful to work 
at a "stretched horizon" = taking e to zero at the end of the computation. Near this stretched 
horizon, one can define a vector orthogonal to the orbits of by 

Vax' = (73) 

where k is the surface gravity at the horizon. Note that 

xVa = --X"X'V„X6 = 0. (74) 

K 
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At the horizon, and become null, and the normalization in Eq. (73) has been chosen so that 
pa ^ Away from the horizon, however, and define two orthogonal directions. 
Then we focus on vector fields of the form 

^« = Rp^ + Tx". (75) 

The corresponding diffeomorphisms are, in a reasonable sense, deformations in the "r-i plane," 
which are known to play a crucial role in the Euclidean approach to black hole thermodynamics. 

Now, given any one-parameter group of diffeomorphisms satisfying boundary conditions, it is 
possible to check that 

{6,6}" = - T2DTi)x" + -^D{TiDT2 - T2DTi)p\ (76) 

This is isomorphic to the standard algebra of diffeomorphisms of the circle or the real line. The 

question is whether the algebra of constraints merely reproduces this DiffS^ or DiffH, algebra, or 
whether it acquires a central extension. Finally, one finds that 



-{TiD^T2 - T2D^Ti) - 2k(TiDT2 - T2DT1) 

K 



(77) 



where terms of order have been omitted. 

This expression has the characteristic three-derivative structure of the central term of a Virasoro 
algebra. We must also compute the surface term ^[{6)6}] of the Hamiltonian to obtain the 
complete expression for the central term in the constraint algebra. This Hamiltonian consists of 
two terms. The first is straightforward to compute, one finds that 

Qai...a„_. = Y^ia,...a„., (2kT - ^D^T^ + O(x'). (78) 

The second term is more complicated, it is shown that it makes no further contribution. Hence we 
obtain a central term 

^[6,6] = {DT,D^T2-DT2D^n) , (79) 

and a centrally extended constraint algebra 

{J[^i],JM}* = + (80) 

A one-parameter group of diffeomorphisms satisfying orthogonality that closes under the brackets 
(76) is then given by 

'^ri = ^ exp I in ^KV + ^ 4(<^(a) - ^ia)v)j | , (81) 

where the ia are arbitrary integers, at least one of which must be nonzero, and the normalization 
has been chosen so that 

{Tm, Tn} = -i{m - n)Tm+n ■ (82) 

Assuming the orthogonality relation, it is easy to see that the algebra (80) is now a conventional 
Virasoro algebra. The microscopic degrees of freedom, whatever their detailed characteristics, must 
transform under a representation of this algebra. But as Strominger observed [13] , this means that 
these degrees of freedom have a conformal field theoretic description, and powerful methods from 
conformal field theory are available to analyze their properties. 
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In particular, we can now use the Cardy formula to count states. If we consider modes of the 
form (81), the central term (79) is easily evaluated: 

iA 

K[Tm, Tn] = --^m'^Sm+nfi, (83) 

where A is the area of the cross section H. The algebra (80) thus becomes 

i{J[Tm],J[Tn]}* = {m- n)J[Tm+n] + -^rrV^Sm+n,o, (84) 



which is the standard form for a Virasoro algebra with central charge 

c A 
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(85) 



The Cardy formula also requires that we know the value of the boundary term J[Tq] of the Hamil- 
tonian. This can be computed from Eq. (78): 

JIT.] = (86) 

The Cardy formula then tells us that for any conformal field theory that provides a representation 
of the Virasoro algebra (84) — modulo certain assumptions — the number of states with a given 
eigenvalue A of J[Tq] grows asymptotically for large A as 



p(A)~exp|2vr^| " ^) } • (8^) 
Inserting Eqs. (85) and (86), we find that 

logp-^, (88) 

giving the expected behavior of the entropy of a black hole. 

Despite the strong desire by Carlip to avoid technical errors this approach has also been criti- 
cized [22] . Another calculation [22] gives zero result for the central charge. The difference is in the 
choice of basis and it is observed that the basis taken by Carlip is singular on the horizon. 



2.5. Carlip III 

A new method was proposed and applied to 2-dimensional dilaton gravity [23]. The main idea 
seems to be that the black hole horizon should contribute directly to the symplectic form and so to 
the Poisson brackets. 

For the discussion of details and criticism we address reader to paper [24]. It seems to be 
appropriate to remind here that a surface contribution to the symplectic form and Poisson brackets 
have been observed and studied in [25, 26]. 



3. Other approaches to statistical derivation of black hole entropy 
3.1. String theory 

String theory becomes related to black holes due to super gravity which is its low energy limit. 
The special interest is attracted to black hole solutions which do not get corrections under quanti- 
zation. These are so-called BPS solitons. String theory methods alow to vary the coupling constant 
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keeping entropy unchanged. For small string coupling constant a black hole behaves as a point 
particle in flat space-time. 

We consider as the first work in this direction the paper by Stroniingcr and Vafa [27], where 
Bckcnstcin-Hawking formula has been derived for one class of 5-dimensional extremal black holes. 
This result has been reached by a calculation of degeneration degree for BPS soliton states. The 
degeneration degree was estimated by topological methods and in leading order the statistical theory 
gave the following formula 



where Qh is the axion charge and Qp — the electrical charge. In the same time proceeding from 
low energy effective action one could derive the Bekenstein-Hawking entropy as 



The obtained results (89) and (90) coincide when charges are large. Qh and Q|^/2 are here integers. 

The essential content of this pioneer work was the exploration of five noncompact dimensions 
and A*" = 4 super symmetry. 5-dimensional space occurred the first space giving a non-trivial result 
for the entropy. To make the horizon area different from zero it was necessary to have nonzero 
values for both charges Qh and Qp- 

Of course, there were a lot of following publications that considerably strengthened first results. 
So, Horowitz and Strominger in article [28] generalized calculations of [27] onto near extremal black 
holes (in the first order deviation from extremal ones) . The excitations of extremal black holes were 
identified with the string states. This allowed to study also a process of black hole evaporation 
in the framework of perturbation theory. It was found that unitarity is preserved in this process. 
5-dimensional black hole was replaced by 6-dimensional black string wrapped around a compactified 
dimension (1-circle). Then the black hole states were described as degrees of freedom inhabiting 
this circle. 

The extremal black hole has 2 sorts of charge: one charge is responsible for the number of degrees 
of freedom and another — for their right momenta. Left momenta are zero in case of the extremal 
black hole. It is supposed that a non-extremal black hole would have also nonzero left momenta. 
The extremal solution with zero right momenta has zero horizon area, i.e it is non-degenerate. But 
if there are nonzero right momenta, then the black string though staying extremal gets a nonzero 
horizon area. 

The dimensional reduction transforms 6-dimensional black string into 5-dimensional black hole 
such as described in paper [27]. If one adds left momenta then the black string becomes not 

extremal and compactifics to a non-extremal black hole. Calculating the number of black string 
states Horowitz and Strominger has obtained a precise correspondence with the black string entropy 
formula. 

The following work by Maldacena and Strominger [29] has allowed to go from 5-dimensional 
extremal black holes to 4-dimensional extremal ones. The methods of works discussed above were 
based on calculations with D-branes. But 4-dimensional black hole with a nonzero horizon area 
could not be constructed of D-branes only. It was necessary to introduce a new object called a 
symmetric 5-brane or Kaluza-Klein monopole and to apply a new technique of calculations. This 
was done in paper [29]. 

Exploiting the standard formula for (1 -|- l)-entropy 




(89) 




(90) 




(91) 
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where Nb (-^f) is a number of boson (fermion) species moving in right direction, is a total 
energy, and L is characteristic length, and also using relations Nb = Np = E = 27rn/L, 

one might get for thermodynamical limit of large n to an independent of L result: 

Sstat = I'K^lQiQ^n . (92) 

The Bekenstein-Hawking entropy is as follows 

Sbh = 2'iT\/Q2QQnm , (93) 

where m is the axion charge of a symmetric 5-brane wrapped around 1" x S"^. As this charge is 
absent in the compactification S^, so Sbh is always proportional to (charge)^ when estimated from 
the leading order low energy effective action. This is quite different from 5-dimensional case when 
it is proportional to (charge)^/'^ . As Sstat is proportional to (charge)'^^'^ , it appears in leading order 
only in 5 dimensions and stays an invisible correction of higher order in 4 dimensions. 

In order to get a nonzero area in 4 dimensions it is required to add m 5-branes, wrapping Yx S^. 
These m 5-branes may be placed anywhere on S^. Each 2-brane has to intersect all m 5-branes 
along S^. A 2-brane can be broken and its ends separate in Y when crossing a 5-brane. Hence the 
Q2 toroidal 2-branes break up into mQ2 cylindrical 2-branes, each of which is bounded by a pair of 
5-branes. The momentum-carrying open strings now carry an extra label describing which pair of 
5-branes they lie in between. The number of species becomes Nb = Np = AmQ2Q&- Inserting this 
into (91) together with E = 2'Kn/L we obtain 

Sstat = 2'K^Q2Q&nm, (94) 

In agreement with the semiclassical result for Sbh- 

For the N = A case there are, in general, 28 electric charges Q and 28 magnetic charges P. 
In this notation 2Q2Q6 = and 2nm = Q^. Duality implies that the entropy depends only on 
P^ , and Q ■ P. The general formula for the Bekenstein-Hawking entropy is 

Sbh = TT^Jp^Q^-iQ-Py. (95) 

For the discussed case the last term vanishes. 
One may ask the following questions: 

• What is the nature of these degrees of freedom responsible for black hole entropy? 

• Where they are located? 

• Is it possible to treat non-extremal black holes by this method? 

They are still open and should be answered if the string theory pretends to be a theory of everything. 

3.2. Loop quantum gravity 

In the Ashtekar approach [30] first a phase space is constructed which describes a space-time 
containing an isolated non-rotating black hole. Then it is canonically quantized. As a result a 
black hole sector in the non-perturbative canonical quantum gravity is obtained. Quantum states 
determining the black hole horizon geometry are separated. Just they are responsible for the 
Bekenstein-Hawking entropy. It is found that these black hole degrees of freedom can be described 
by Chern-Simons field theory on the horizon. It is shown that the statistical mechanical entropy is 
proportional to the horizon area for large non-rotating black holes. The constant of proportionality 
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depends upon the Immirzi parameter [31], which fixes the spectrum of the area operator in loop 
quantum gravity; an appropriate choice of this parameter gives the Bekenstein-Hawking formula 

S = yl/4£pp The same choice of the Immirzi parameter also gives the similar result for black holes 
carrying electric or dilatonic charge. An important advantage of this method is that black holes 
considered are not necessary extremal or near extremal, and also they are certainly 4-dimensional 
(not 5-dimensional !). 

Non-perturbative techniques of the Ashtekar approach have led to a quantum theory of geometry 
in which operators corresponding to lengths, areas and volumes have discrete spectra [32]. Of 
particular interest are the spin network states associated with graphs in 3-space with edges labelled 
by spins j = ^, 1, . . . and vertices labelled by intertwining operators. If a single edge punctures a 
2-surface transversely, it contributes an area proportional to 

A sector of the classical phase space corresponding to an isolated, uncharged, non-rotating black 
hole can be described starting from a manifold with a boundary. This space-time manifold has the 
outer boundary I and the inner boundary asH. I corresponds to the infinitely far region and H — 
to the black hole horizon. Hamiltonian variables are a soldering form for SL(2, C) spinors 

and an SL(2, C) connection AaA^ ■ Metric Qah = o"^"^ (^hAA' is the Lorentzian space-time metric and 
AaA^ is the self-dual connection operating on unprimed spinors only. On I fields are required to 
satisfy the standard asymptotically flat boundary conditions. The conditions on H, on the other 
hand, are specially elaborated by Ashtekar's group. The key requirements are: i) H should be a null 
surface with respect to the metric gab', ii) On a 'finite patch' A of TC, the area of any cross-section 
should be a constant. As, the Weyl spinor be of Petrov type 2-2 and its only non-zero component, 
^'2, should be given by *2 = "^t^/As; and, hi) the 2-fiats on A, orthogonal to the two principal 
null directions of the Weyl tensor span 2-spheres and the pull-back of the connection Aa to these 
2-spheres should be real. 

In short, condition ii) implies that there is no gravitational radiation falling into A ('isolated' 
black hole) while the first part of iii) implies that it is 'non-rotating'. These boundary conditions 
have been taken from the Schwarzschild black hole horizon. But here space-time is not required to 
be static and gravitational waves are admitted far from the black hole. Therefore, the phase space 
is infinite dimensional. Nevertheless, these boundary conditions are sufficient for the existence of 
a surface integral providing differentiability of the action, so the Regge-Teitelboim criterion [16] is 
fulfilled. The required surface term is just the Chern-Simons action: 

where ^^i^ = 2a[a'^^ ^b]A'^ , ^'^^ FabA^ is the curvature of the connection A. In the Hamiltonian 
formalism following from (96) canonical variables are the restrictions of S and A to the spatial 
hypersurface M with a boundary S. Here a difficulty arises due to the fact that the restriction to 
M of the self-dual connection ^ is a complex valued SU(2) connection. Then it is better to make a 
transition to real variables: Aa = Ta — iKa, where T is the 3-dimensional spin connection compatible 
with the triad field and K is the extrinsic curvature of M. Here the Immirzi [31] parameter 7 appears. 
It is a positive real number and real phase space variables are lAa := Fa — jKa and '^ab '■= {^/l)^ab- 
Then the boundary conditions imply that the phase space consists of real fields (lAa,'^ab) that are 
asymptotically flat at inflnity and satisfy the following condition at S: 

^zr = -^^4^ (97) 

where underlining means pull-backs to S. The restriction of lAa to S yields a reducible connection 
satisfying DaV = for some radial internal vector r. It is suitable to fix r on 5 using the SU(2) 
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gauge freedom. Then the gauge group on the boundary is reduced to U(l) and only the r component 
of (97) is non-trivial. 

On the real phase space, the symplectic structure corresponding to (96), in addition to the 
familiar volume term, contains a surface contribution (sec [25] for the first evidence of a surface 
term in Ashtekar's symplectic form) which coincides with the symplectic structure of the Chern- 
Simons theory: 

n = / Tr\6"^A 5'^.' - 6"^' A S'^A] - —^f-^ I TV \5^A A 6^A'] . (98) 
SttGJm Wtt^jG Js 

There are bulk degrees of freedom, corresponding to gravitational waves far away from A, which 
should not be taken into account as genuine black hole degrees of freedom. The boundary degrees 
of freedom describing the geometry of the horizon S are degrees of freedom 'living on the horizon' 
and responsible for the entropy. 

In the classical theory the bulk and boundary degrees of freedom cannot be separated: all fields 
on S are determined by fields in the interior of M as they are required to be continuous. However, 
in the quantum theory, the fields describing geometry become discontinuous, and the fields on S are 
no longer determined by fields in M; in this case there are independent degrees of freedom 'living' 
on the boundary. These surface degrees of freedom are the ones that account for black hole entropy 
in this approach. 

To quantize the theory a Hilbert space of 'volume' states and a Hilbcrt space of 'surface' 
states are constructed, and then constraints are imposed on Ti^ ®1-L^ to obtain the space of physical 
states. HX is taken to consist of square-integrable functions on the space of generalized SU(2) 
connections on M modulo gauge transformations that are identity on S. The form of the Hilbert 
space of surface states is motivated by the fact that in the quantum theory it is required to 
impose the boundary condition (97) as an operator equation. That is, given a spin network state 
'^v in and a state s in ■> the quantum version of the r component of equation (97) should 
read 

(l(8)^Zaf, -r + t^,- r(g)l) ^v®'i>s = 0. (99) 

The structure of this equation implies that ^'y and ^'g should be eigenstates of • r and • r 

respectively. The 'polymer nature' of quantum geometry in M implies that eigenvalues of • r 
are distributional, given by 

87r4i jiS\x,pi)vab ill (100) 

i 

for some points pi on S, where ji are half-integers, 6^ is the delta distribution on S, rjab the Levi- 
Civita density on S and £pi the Planck length. Therefore, (99) implies that the surface states *5 
have support only on generalized connections that are everywhere flat except at a finite number of 
points Pi. It turns out that such generalized connections can be identified with ordinary connections 
with distributional curvature. Since the surface symplectic structure is that of Chern-Simons theory, 
for any fixed choice 

'P = {{Pl,jpi),---,{Pn,jpJ} 

of points in S labelled by spins, should have a subspace given by the space of states of U(l) 
Chern-Simons theory on a sphere with punctures p labelled by spins jp. The total space is the 
direct sum of these subspaces. 

Note now that (k/2Tr)F_ is the generator of internal rotations in Chern-Simons theory. Thus, 
the meaning of (99) turns out to be rather simple: it ensures that the volume and surface states are 
'coupled' in precisely the correct way so that the total state is invariant under U(l) internal rotations 
at S. The remaining constraints require that the states be invariant under diffeomorphisms of M 
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that leave S invariant and under motions generated by the Hamiltonian constraint smeared with 
any lapse field that vanishes at S. Thus for each set V of finitely many punctures p labelled by 
spins jp, there is a subspace of volume states having a basis given by open spin networks whose 
edges intersect S only at these punctures, where they are labelled by the spins jp. Similarly there is 
a subspace Hj, consisting of quantum states of U(l) Chern-Simons theory on the punctured surface 
S. The total physical Hilbert space is given by: 

tphy — 



Gauge 

where 'Gauge' means internal SU(2) rotations that reduce to U(l) on S, diffeomorphisms preserving 
S, and the motions generated by the Hamiltonian constraint. To proceed further it is necessary 
to make an assumption about the quantum dynamics: let there is at least one solution of this 
constraint in Tip Ti-j, for any set V of punctures labelled by spins. 

It is required to consider only states of the horizon of a black hole with area A. Thus the 'volume' 
states are traced over to construct a density matrix pbh describing a maximal-entropy mixture of 
surface states for which the area of the horizon lies in the range As^l'^i- The statistical mechanical 
black hole entropy is then given by 5'bh = — Trpbh Inpbh- This can be computed by counting states: 
^bh = In A^bh where A^bh is the number of Chcrn- Simons surface states satisfying the area constraint. 

The eigenvalues of the area operator are given by the following formula [32]: 

87r74iEVS(^+l) > (101) 
p 

where jp are the spins labelling the punctures. Using this and the result from Chern-Simons theory 
that for a large number of punctures the dimension of TCj, grows as 

dimHl ~ n + 1) > (102) 

it is possible to calculate the entropy. For large A it is given by 

e 70 . In2 
'^hh = TTo As, 7o 



Thus, in the limit of large area, the entropy is proportional to the area of the horizon. If 7 is put 
equal to 70, then the calculation based on statistical mechanics gives a result coinciding with the 
Bekenstein-Hawking formula. 

As it is stressed by authors [30], the advantages of their method are 1) that black holes considered 
are not necessary extremal or near extremal, and 2) that dimension of space-time has just the 
observable value. Nevertheless, as one can see from above, Ashtekar's approach to quantum gravity 
faces with some principal difficulties. First, they seem to be coming from the complex structure of 
the phase space, second, from the fact that space-time dimension here can not be easily changed in 
Ashtekar's formalism. Probably, other difficulties, such as the absence of proof that the Hamiltonian 
constraint solutions really exist, and the presence of a free parameter (Immirzi parameter 7), are 
consequences of these two. 
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4. Conclusions 

• There is a strong belief that the problem of statistical derivation of the Bekenstein-Hawking 
formula for black hole entropy can be solved at the present level of our physical knowledge. 

• It is very probable that the path to solution consists in exploiting conformal invariance near 
the black hole horizon. 

• But now this problem is far from being completely solved. 

• Every successful approach to the problem could help us to understand what quantum gravi- 
ty is. 

It would be unrealistic to try to cover all approaches to statistical derivation of the black hole 

entropy here, as an example let us mention rather radical ideas appeared quite recently in works 
by V.V. Kiselev [33]. We would like to address interested readers to more detailed reviews of the 
subject [12, 34, 35]. 
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